Optimal cloning of coherent states with a linear amplifier and beam splitters 
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A transformation achieving the optimal symmetric N — > M cloning of coherent states is presented. 
Its implementation only requires a phase-insensitive linear amplifier and a network of beam splitters. 
An experimental demonstration of this continuous-variable doner should therefore be in the scope 
of current technology. The link between optimal quantum cloning and optimal amplification of 
quantum states is also pointed out. 

PACS numbers: 03.67.-a, 03.65. Bz, 42.50.Dv, 89.70.+C 



Quantum systems cannot be cloned exactly pi], but 
only approximately. Finding the optimal approximate 
quantum cloning transformation has been a fundamen- 
tal issue in quantum information theory for the last five 
years. In quantum cryptography, for instance, this prob- 
lem happens to be strongly related to the assessment 
of security . Cloning has been extensively studied to 
date for discrete quantum variables, such as quantum bits 
|^-^| or ci-level systems [fl0|-[i"2|, since quantum informa- 
tion theory was initially developed for these kinds of sys- 
tems. Recent progress has shown, however, that contin- 
uous spectrum systems might be experimentally simpler 
to manipulate than their discrete counterparts in order 
to process quantum information (see, e.g., |l3|Jl^ ]). 

Stimulated by this progress, we investigate in this Let- 
ter the possibility of achieving the cloning of continuous- 
variables quantum information. Commonly, a distinc- 
tion is made between universal cloning, if the set of in- 
put states contains all possible states for a given Hilbert 
space dimension, and state- dependent cloning, if the in- 
put states are restricted to a certain set which does not 
contain all possible states. For any Hilbert space dimen- 
sion, the optimal universal doner fl(i|-|r^| that clones all 
possible input states equally well can be constructed from 
a single family of quantum circuits fl5f| . This univer- 
sal doner reduces to a classical probability distributor in 
the continuous limit. Besides the universal doner, quan- 
tum cloning of continuous- variable systems has been con- 
sidered first in a state-dependent context. In Ref. [ jl6| , 
the duplication of coherent states was investigated, and 
an explicit transformation that is covariant under dis- 
placement and rotation in phase space was derived. This 
transformation therefore clones all coherent states with 
the same fidelity (F = 2/3), although it is not univer- 
sal, strictly speaking, as its cloning fidelity is lower for 
other classes of states such as squeezed states. The op- 
timality of this continuous-variable cloning transforma- 
tion was then proven in Ref. [l7| . More generally, it was 
shown that, if one attempts to produce M clones from 
N original replicas of a coherent state \a) (M > N) with 



an equal fidelity for all ev's, the so-called N-to-M cloning 
transformation must result in an additional noise on both 
quadratures of each of the M outputs which has a mini- 
mum variance 
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where the vacuum noise on a quadrature is denoted as 
Ax 2 ac = 1/2 (H = 1). The corresponding maximum 
N -> M cloning fidelity is F NM = M N/(MN + M — N). 
However, finding the optimal N — > M cloning transfor- 
mation and proving that it actually achieves this maxi- 
mum fidelity was an open problem. 

The present Letter resolves this question. We use 
the Heisenberg picture in order to derive explicitly an 
N — ► M symmetric cloning transformation that at- 
tains Eq. ([!]). Remarkably, it appears that implement- 
ing this transformation only requires a phase-insensitive 
linear amplifier and a set of beam splitters. Let = 



a 



® |0) 



®M — N 



® |0) z denote the initial joint state of 



the N input modes to be cloned (prepared in the coher- 
ent state |a)), the additional M — N blank modes, and 
an ancillary mode z. The blank modes and the ancilla 
arc assumed to be initially in the vacuum state |0). Let 
{xk,Pk} denote the pair of quadrature operators associ- 
ated with each mode k involved in the cloning transfor- 
mation: k — ... N — 1 refers to the N original input 
modes, and k = N . . . M — 1 refers to the additional 
blank modes. Cloning can be thought of as some unitary 
transformation U acting on |^), and resulting in a state 
\*ff") = U\if>) such that the M modes are left in the same 
(mixed) state which is maximally close to |or). Alterna- 
tively, in the Heisenberg picture, this transformation can 
be described by a canonical transformation acting on the 
operators {x k ,p k }: 



x k U, 
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Pk u, 



(2) 



while leaving the state |^) invariant. We will work in the 
Heisenberg picture and use the above notation through- 



1 



out this paper, with x'£ denoting the clones (i.e., the out- 
put modes of the cloning circuit except the ancilla z), 
because cloning turns out to be much simpler to describe 
from that point of view. We will now impose several re- 
quirements on transformation (Q) that translate the ex- 
pected properties for an optimal cloning transformation. 
First, we require that the M output modes have the de- 
sired mean values: 



[Xk 



= [a\x \a) 



(Pk") = (a\p \a), 



(3) 



for k = ... M — 1. Roughly speaking, this means that 
the state of the clones is centered on the original coher- 
ent state. Our second requirement is covariance with 
respect to rotation in phase space. Coherent states have 
the property that the quadrature variances are left in- 
variant by complex rotations in phase space. So, for any 
input mode k of the cloning transformation and for any 
operator Vk — cxk + dpk (where c and d are complex 
numbers satisfying |c| 2 + \d\ 2 = 1), the error variance a% k 
is the same: 



(4) 



We impose this property to be conserved through the 
cloning process. Taking optimality into account, Eq. (Q), 
rotational covariance yields: 
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where 



dp'l- Our third requirement is, of 



course, the unitarity of the cloning transformation Q). 
In the Heisenberg picture, this is equivalent to demand- 
ing that the commutation relations are preserved through 
the evolution: 



,Xk"] = W',Pk"] = 0, 
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for j, k = . . . M — 1 and for the ancilla. 

Let us first focus on the continuous-variable duplica- 
tion (N — 1, M — 2). A simple transformation obeying 
the three conditions mentioned above is given by: 



xo _ xo + — + — 



, Pi Pz 
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X\ 



V2 V2' 
x' z = x + v2x z , 



Pi = Po 



Pz = -Po 



Pi Pz 



(7) 



This transformation clearly conserves the commutation 
rules, and yields the expected mean values (xq,pq) for 
the two clones (modes 0" and 1"). Also, one can check 
that the quadrature variances of both clones are equal 
to 2Ax^ ac , in accordance with Eq. (||). This transfor- 
mation actually coincides with the Gaussian cloning ma- 
chine introduced in Ref. Jl6| . Interestingly, we note here 
that the state in which the ancilla z is left after cloning 
is centered on (xq,— Pa), that is the phase- conjugated 



state | a*). This means that, in analogy with the uni- 
versal qubit cloning machine ||], the continuous- variable 
doner generates an "anticlone" (or time-reversed state) 
together with the two clones. 

Now, let us show how this duplicator can be imple- 
mented in practice. Equation (Q) can be interpreted as a 
sequence of two canonical transformations: 
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where = (xk + ipk)/V% and 
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note the annihilation and creation operators for mode fc. 
As shown in Fig. |], the interpretation of this transforma- 
tion becomes then straightforward: the first step (which 
transforms ao and a z into a' and a' z ) is a phase-insensitive 
amplifier whose (power) gain G is equal to 2, while the 
second step (which transforms a' and a\ into Oq and a'{) 
is a phase-free 50:50 beam splitter fTq ]. Clearly, rota- 
tional covariance is guaranteed here by the use of a phase- 
insensitive amplifier . As discussed in Ref. [fl9t , the ancilla 
z involved in linear amplification can always be chosen 
such that (a z ) = 0, so that we have (ad) = (a'{) = (a ) 
as required. Finally, the optimality of our doner can be 
confirmed from known results on linear amplifiers. For 
an amplifier of (power) gain G, the quadrature variances 
of a z are bounded by |19|: 



> (G-l)/2. 



(9) 



Hence, the optimal amplifier of gain G = 2 yields a\ — 
1/2, so that our cloning transformation is optimal. 
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FIG. 1. Implementation of a 1 — > 2 continuous- variable 
cloning machine. LA stands for linear amplifier, and BS rep- 
resents a phase-free 50:50 beam splitter. 

Let us now consider the N — > M cloning transfor- 
mation. In order to achieve cloning, energy has to be 
brought to the M — N blank modes in order to drive 
them from the vacuum state into a state which has the 
desired mean value. We will again achieve this operation 
with the help of a linear amplifier. From Eq. (^), we see 
that the cloning induced noise essentially originates from 
the amplification process, and grows with the amplifier 
gain. So, we should preferably amplify as little as possi- 
ble. Loosely speaking, the cloning procedure should then 
be as follows: (i) symmetrically amplifying the N input 
modes by concentrating them into one single mode, which 
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is then amplified; (ii) symmetrically distributing the out- 
put of this amplifier among the M output modes. As 
we will see, a convenient way to achieve these concentra- 
tion and distribution processes is provided by the Dis- 
crete Fourier Transform (DFT). Cloning is then achieved 
by the following three-step procedure (see Fig. First 
step: a DFT (acting on N modes), 
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exp(ikl2ir / N) ai, 



(10) 



with fc = . . . N — 1. This operation concentrates the en- 
ergy of the N input modes into one single mode (renamed 
ao) and leaves the remaining N — 1 modes (a[ . . . a / N _ 1 ) 
in the vacuum state. Second step: the mode ao is am- 
plified with a linear amplifier of gain G = M/N . This 
results in 



(11) 



Third step: amplitude distribution by performing a DFT 
(acting on M modes) between the mode a' and M — 1 
modes in the vacuum state: 




^ Af-l 
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exp(ikl2ir/M) a[, 



(12) 



with k = . . . M - 1, and a- = <Zj for i = N . . . M - 1. 
The DFT now distributes the energy contained in the 
output of the amplifier among the M output clones. 
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FIG. 2. Implementation of an N — > M continuous-variable 
cloning machine. C stands for the amplitude concentration 
operation, while D refers to the amplitude distribution. Both 
can be achieved by using a DFT, or, alternatively, an inverse 
"TV-splitter" and an "M-splitter" (in which case we shift the 
indices by one in the text, i.e., k = 1 . . . M). 

It is readily checked that this procedure meets our 
three requirements, and is optimal provided that the am- 
plifier is optimal, that is a\ t = [(M/N) - l]/2. The 
quadrature variances of the M output modes coincide 
with Eq. (||). As in the case of duplication, the quality 



of cloning decreases as u 2 a increases, that is cloning and 
amplifying coherent states are two equivalent problems. 
For 1 — ► 2 cloning, we have seen that the final amplitude 
distribution among the output clones is achieved with a 
single beam splitter. In fact, any unitary matrix such 
as the DFT used here can be realized with a sequence 
of beam splitters (and phase shifters) [^0| . This means 
that the N — » M cloning transformation can be imple- 
mented using only passive elements except for a single 
linear amplifier. 

We will now explicitly give the simplest beam splitter 
combination that enables the above transformation. For 
convenience, let us now use the indices k = 1 . . . N for the 
N original input modes a^, and k = N + 1 . . . M for the 
additional blank modes at- With an ideal (phase- free) 
beam splitter operation acting on two modes Ck and q, 
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(13) 



we define a matrix Bki{9) which is an M-dimensional 
identity matrix with the entries Iki, Iik, and la re- 
placed by the corresponding entries of the above beam 
splitter matrix. Now we can define a sequence of beam 
splitters acting on M modes ("M-splitter" pi) as 



U(M) = B M -i 



M (sin- 1 -L ) B 
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The individual beam splitters in Eq. ( |14|) depend only on 
their reflectivity/transmittance parameter 9. In order to 
concentrate the N identical inputs, we send them now 
through an inverse TV-splitter, 
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■ U\N) 



ai a 2 



a N ) ■ (15) 



Again, we end up with one mode (renamed ai) hav- 
ing non-zero mean value and N — 1 modes (a' 2 ■ ■ ■ a' N ) 



in the vacuum state. 

y/M/Ntt! 



After amplifying mode ai, 



a. 



ii -t- -J 'M/N — 1 a\ , etc., a final M-splitter oper- 
ation yields the output clones: 



(a'l a!{ ■■■ a%) =U(M)(a' 1 a' 2 
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(16) 



with = ^ for i = N + 1 . . . M. 

Since the amplification produces extra noise, our 
cloning circuits used as little amplification as possible. 
However, rather surprisingly, by first amplifying each in- 
put copy k = 1 . . . N individually, 



(17) 



a circuit can also be constructed that yields optimum 
fidelities. In the next step, the amplified modes are 
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each sent together with M — 1 vacuum modes bfci, 
bk,2,---,bk,M-i through an M -splitter 
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The NM output modes after this operation can be writ- 
ten as 
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where / = 1 . . . M. The noise in each M-splitter output 
coming from the M — 1 vacuum inputs is represented by 
mode dk,i having zero mean value and quadrature vari- 
ances of (M — 1)/2M. The final step now consists of M 
inverse ./V-splitters acting on all modes with the same in- 
dex I, i.e., the ./V modes a', «, and the N modes a', 2 , etc. 
The output modes at each TV-splitter, 



ei.i 



CN-1,1 ) = 



contain only noise except for one mode, 
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Again, all M clones are optimal, although additional 
noise has been introduced at the intermediate steps which 
results in M(N — 1) "waste" output modes. However, 
this particular circuit points out that N — > M cloning of 
coherent states is effectively a "classical plumbing" pro- 
cedure distributing classical amplitudes. 

Finally, we note that for squeezed-state inputs rather 
than coherent states, the transformations and circuits 
presented require all auxiliary vacuum modes (the blank 
modes and the ancillary mode z) be correspondingly 
squeezed in order to maintain optimum cloning fideli- 
ties. This means, in particular, that the amplifier mode 
z needs to be controlled which requires a device differ- 
ent from a simple phase-insensitive amplifier, namely 
a two-mode parametric amplifier. One can say that 
the cloning machine capable of optimum cloning of all 
squeezed states with fixed and known squeezing then op- 
erates in a non-universal fashion with respect to all pos- 
sible squeezed states at the input |l6| , |i"7| ] . 

In summary, an optimal N-to-M continuous-variable 
cloning transformation for coherent states has been 
derived, which attains the maximum cloning fidelity 
F NjM = MN/(MN + M - N). A possible experimen- 
tal implementation of this doner has been proposed. We 
trust that this implementation should be achievable with 
current technology since it only requires a single linear 
amplifier and N + M — 2 beam splitters. In Ref. pj| , 
an alternative one-to-two cloning scheme has been pro- 
posed based on three non-degenerate optical parametric 
amplifiers, and its experimental realization is currently 



underway at Roma University. Finally, we pointed out 
the link between the quality of the best doner and the 
minimum noise induced by the amplification of a quan- 
tum state, emphasizing that spontaneous emission is here 
again the mechanism that prevents the perfect cloning of 
quantum states of light . 
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